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We formulate a universal characterization of the many-particle quantum entanglement in the 
ground state of a topologically ordered two-dimensional medium with a mass gap. We consider a 
disk in the plane, with a smooth boundary of length L, large compared to the correlation length. 
In the ground state, by tracing out all degrees of freedom in the exterior of the disk, we obtain a 
marginal density operator p for the degrees of freedom in the interior. The von Neumann entropy 
S[p) of this density operator, a measure of the entanglement of the interior and exterior variables, 
has the form S[p) — aL — 7 + - • •, where the ellipsis represents terms that vanish in the limit L — > 00. 
The coefficient a, arising from short wavelength modes localized near the boundary, is nonuniversal 
and ultraviolet divergent, but —7 is a universal additive constant characterizing a global feature of 
the entanglement in the ground state. Using topological quantum field theory methods, we derive 
a formula for 7 in terms of properties of the superselection sectors of the medium. 
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In a quantum many-body system at zero temperature, 
a quantum phase transition may occur as a parameter 
varies in the Hamiltonian of the system. The two phases 
on either side of a quantum critical point may be charac- 
terized by different types of quantum order, the quantum 
correlations among the microscopic degrees of freedom 
have qualitatively different properties in the two phases. 
Yet in some cases, the phases cannot be distinguished by 
any local order parameter. 

For example, in two spatial dimensions a system with a 
mass gap can exhibit topological order [1]. The quantum 
entanglement in the ground state of a topologically or- 
dered medium has global properties with remarkable con- 
sequences. For one thing, the quasiparticle excitations of 
the system (anyons) exhibit an exotic variant of indistin- 
guishable particle statistics. Furthermore, in the infinite- 
volume limit the ground-state degeneracy depends on the 
genus (number of handles) of the closed surface on which 
the system resides. 

While it is clear that these unusual properties emerge 
because the ground state is profoundly entangled, up un- 
til now no firm connection has been established between 
topological order and any quantitative measure of en- 
tanglement. In this paper we provide such a connection 
by relating topological order to von Neumann entropy, 
which quantifies the entanglement of a bipartite pure 
state. 

Specifically, we consider a disk in the plane, with a 
smooth boundary of length L, large compared to the 
correlation length. In the ground state, by tracing out 
all degrees of freedom in the exterior of the disk, we 
obtain a marginal density operator p for the degrees 
of freedom in the interior. The von Neumann entropy 
S{p) = — trplogp of this density operator, a measure of 
the entanglement of the interior and exterior variables, 
has the form 



where the ellipsis represents terms that vanish in the limit 
L — > 00. The coefficient a, arising from short wavelength 
modes localized near the boundary, is nonuniversal and 
ultraviolet divergent [2] , but —7 (where 7 is nonnegative) 
is a universal additive constant characterizing a global 
feature of the entanglement in the ground state. We call 
—7 the topological entanglement entropy. 

This universal quantity reflects topological properties 
of the entanglement that survive at arbitrarily long dis- 
tances, and therefore can be studied using an effective 
field theory that captures the far-infrared behavior of 
the medium, namely a topological quantum field theory 
(TQFT) that describes the long-range Aharonov-Bohm 
interactions of the medium's massive quasiparticle exci- 
tations. We find 



7 = logV , 



(2) 



where I? > 1 is the total quantum dimension of the 
medium, given by 



(3) 



S{p) — aL — + 



(1) 



here the sum is over all the superselection sectors of the 
medium, and da is the quantum dimension of a particle 
with charge a. 

Any abelian anyon has quantum dimension d = 1; 
therefore, for a model of abelian anyons, is simply 
the number of superselection sectors. Thus for a Laugh- 
lin state [3] realized in a fractional quantum Hall system 
with filling factor u = 1/q where q is an odd integer, we 
have V = ^/q. For the toric code [4] , which has four sec- 
tors, the topological entropy is 7 = log 2, as has already 
been noted in [5]. 

However, nonabelian anyons have quantum dimension 
greater than one. The significance of da (which need not 
be a rational number) is that the dimension Naaa---a of 
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the fusion vector space spanned by all the distinguishable 
ways in which n anyons of type a can be glued together 
to yield a trivial total charge grows asymptotically like 
the nth power of da- For example, in the SU{2)k Chern- 
Simons theory, we have 




(a) (b) 

FIG. 1. (a) The plane is divided into four regions, labeled 
A, B, C, D, that meet at double and triple intersections, (b) 
Moving the triple intersection where B,C',D meet deforms 
the regions as shown. 

To justify computing the entropy using effective field 
theory, we require that the boundary have no "sharp" 
features that might be sensitive to short-distance physics; 
yet for a boundary drawn on a lattice, sharp corners are 
unavoidable. Furthermore there is an inherent ambigu- 
ity in separating the term that scales with the length L 
from the constant term. We can circumvent these diffi- 
culties by the following construction. We now divide the 
plane into four regions, all large compared to the corre- 
lation length, labeled A, B, C, D as in Fig. la. Let Sa 
denote the von Neumann entropy of the density operator 
PA that is obtained from the ground state by tracing out 
the degrees of freedom outside region A, let Sab denote 
the von Neumann entropy of the density operator pab 
obtained by tracing out the degrees of freedom outside 
region AB = AU B, etc. Then we define the topological 
entropy S'topo as 

5'topo = Sa + Sb + Sc — Sab — Sbc — Sac + Sabc ■ 

(5) 

This linear combination of entropies has been strategi- 
cally chosen to ensure that the dependence on the length 
of the boundaries of the regions cancels out. For exam- 
ple, the term proportional to the length of the double 
intersection of A and D appears in Sa and Sabc with 
a -I- sign, and in —Sab and —Sac with a minus sign. 
Similarly, the double intersection of A and B appears in 
Sa and Sb with a -I- sign, and in —Sac and in —Sbc 
with a minus sign. (The observation that the ultraviolet 
divergent terms cancel in a suitably constructed linear 
combination of entropies has also been exploited in [6] 
and applied there to (l-l-l)-dimensional systems.) 



Assuming the behavior eq. (1) in each term, we find 
S'topo — ~7- But the advantage of defining iStopo using a 
division into four regions is that we can then argue per- 
suasively that S'topo is a topological invariant (dependent 
only on the topology of how the regions join and not on 
their geometry) and a universal quantity (unchanged by 
smooth deformations of the Hamiltonian unless a quan- 
tum critical point is encountered). 

To see that S'topo is topologically invariant, first con- 
sider deforming the boundary between two regions, far 
from any triple point where three regions meet. Deform- 
ing the boundary between C and D, say, has no eflFect 
on regions A, B, and AB] therefore if all regions are 
large compared to the correlation length, we expect the 
changes in Sa, Sb, and Sab to all be negligible. Thus 
the change in S'topo can be expressed as 

AStopo = (ASabc - ASbc) - (ASac - ASc) . (6) 

We expect, though, that if the regions are large compared 
to the correlation length, then appending region A to BC 
should have a negligible effect on the change in the en- 
tropy, since A is far away from where the deformation is 
occuring; similarly, appending ^ to C should not affect 
the change in the entropy. Thus both terms on the right- 
hand side of eq. (6) vanish, and Stopo is unchanged. The 
same reasoning applies to the deformation of any other 
boundary between two regions. 

Next consider deforming the position of a triple point, 
such as the point where B, C, and D meet as in Fig. lb. 
Again we may argue that Sa is unchanged by the defor- 
mation. We recall that for a bipartite pure state (like the 
ground state), the marginal density operators for both 
subsystems have the same nonzero eigenvalues and there- 
fore the same entropy; thus Sabc = Sd and Sbc = Sad- 
We see that the change in Stopo can be expressed as 

AStopo = (ASb - ASab) 
+(ASc - ASac) 

+ i^SD-ASAD) ■ (7) 

All three terms on the right-hand side of eq. (7) van- 
ish because appending A does not aflFect the change in 
the entropy. The same reasoning applies when any other 
triple point moves; we conclude that Stopo is unchanged 
by any deformation of the geometry of the regions that 
preserves their topology, as long as all regions remain 
large compared to the correlation length. 

Now, what happens to Stopo as the Hamiltonian of 
the system is smoothly deformed? We assume that the 
Hamiltonian is a sum of local terms, and that the corre- 
lation length remains finite during the deformation (and 
in fact that the correlation length stays small compared 
to the size of regions A, B,C, D). If the Hamiltonian 
changes locally in a region far from any boundary, then 
this change has a negligible effect on the ground state 
in the vicinity of the boundary, and therefore does not 
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affect 5'topo- If the Hamiltonian changes locally close to 
a boundary, we can exploit the topological invariance of 
S'topo to first move the boundary far way, then deform 
the Hamiltonian, and finally return the boundary to its 
original location. Thus wc see that S'topo is a universal 
quantity characteristic of a particular kind of topologi- 
cal order, which remains invariant if no quantum critical 
point is encountered as the Hamiltonian varies. 




FIG. 2. The planar medium is glued at spatial infinity to 
its time-reversal conjugate, and wormholes are attached that 
connect the two conjugate media at the locations of the triple 
intersections, creating a sphere with four handles. Each re- 
gion, together with its image, becomes a sphere with three 
punctures, and each union of two regions, together with its 
image, becomes a sphere with four punctures. The punctures 
carry charges labeled a, b, c, d. Anyons that wind around a 
cycle enclosing a wormhole throat as shown detect a trivial 
charge. 

To facilitate the computation of S'topo; it is convenient 
to imagine joining together the planar medium we wish 
to study with its time-reversal conjugate. We glue to- 
gether the medium and its conjugate at spatial infinity, 
and then attach "wormholes" that connect the two planes 
at the positions of the four triple intersections, as in- 
dicated in Fig. 2. The resulting closed surface has the 
topology of a sphere with four handles. If an isolated 
planar medium is punctured, then massless chiral modes 
propagate around the edge of the puncture, but the edge 
states of the medium and its conjugate have opposite 
chirality, so the edge states acquire masses when the two 
surfaces are coupled; therefore, the wormholes can be cre- 
ated adiabatically without destroying the mass gap. No 
anyons are produced during this adiabatic process, so 
that the mouth of each wormhole carries trivial anyonic 
charge. 

The boundaries that separate regions in the plane and 
in its double can be joined through the wormholes as 
in Fig. 2; then each region of the doubled surface has 
the topology of a sphere with three punctures, and each 
union of two adjacent regions becomes a sphere with four 
punctures. The topological entanglement entropy of the 
medium and its conjugate are both equal to S'topo, so that 



the topological entanglement entropy of the doubled sur- 
face is twice Stopo- The entropy of a region depends only 
on its topology, so for the doubled surface we have 

2Stopo — 4S3 — 3S4 , (8) 

where S3 denotes the entropy for the sphere with three 
punctures and S4 denotes the entropy for the sphere with 
four punctures. 

The quantities S3 and S4 can be computed using the 
appropriate effective field theory, a TQFT [7]. We use 
the property that no charge is detected by an anyon that 
winds around the throat of a wormhole. A cycle that 
encloses a puncture in the double of (say) region A is 
complementary to a cycle that winds around the worm- 
hole throat; it follows that the puncture carries charge a 
with probability 

Pa^\S^\'=dl/V\ (9) 

where is the topological S-matrix of the TQFT, and 
1 denotes the trivial charge. To find the joint probabil- 
ity distribution pabc governing the charges a, 6, c on the 
punctures of the sphere with three punctures, we may 
use standard TQFT methods to compute the probability 
Pab-^c that when charges a and b fuse the total charge is 
c. The result is 

Pab-^c = Nabcdc/dadb , (10) 

where Nabc is the dimension of the fusion vector space 
spanned by all the distinguishable ways in which charges 
a, b, and c can fuse to yield trivial total charge; it follows 
that 

Pabc = PaPb ■ Pab^c = Nabcdadbdc/T)'^ ■ (H) 

Evaluating the entropy in the basis in which each punc- 
ture has a definite charge, and summing over all the dis- 
tinguishable fusion states that occur for specified values 
of the charges, we find 

r. Pabc , / Pabc \ 

a,o,c fi—l 

= 4l0gP- ^Pabc'^Og{dadbdc) 
a,b,c 

= UogV- 3 J^Palog da . (12) 

a 

For the sphere with four punctures, a similar calculation 
yields 

Pabcd = PaPbPc ■ Pabc^d = N abcddadbdcdd/V^ (13) 

and 

S4 = 61ogl?-4^p,logd, . (14) 

a 

Plugging into eq. (8), we find 
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^topo = 253 - -54 = - log I? = -7 . (15) 

Eq. (15) is our main result. Note that it follows if we 
use eq. (1) to evaluate the entropy of each region, since 
7 appears four times in the expression for 5topo with a 
negative sign and three times with a positive sign. We 
also observe that 5topo actually depends on the topol- 
ogy of the regions A, B, C. For example, consider the 
arrangement shown in Fig. 3, in which B and AC both 
have two connected components, and ABC is not sim- 
ply connected. Since regions B, AC, and ABC each 
have boundaries with two components, now 7 appears 
six times with a negative sign and four times with a pos- 
itive sign, so that 5topo = — 27. 
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FIG. 3. If the regions A,B,C have the topology shown, 
then Stopo = —27. 

Using a different approach, we can formulate a sim- 
pler but more heuristic derivation of the formula for 7. 
First we write the marginal density operator p for the 
disk as p = e~^^ . This is just a definition of H and 
has no other content in itself; furthermore the param- 
eter T = j3~^ \s arbitrary — so we are free to choose 
it to be small compared to the bulk energy gap of the 
two-dimensional medium. Now we make a natural but 
nontrivial assumption: that H can be regarded as the 
Hamiltonian of a (1-1- l)-dimcnsional conformal field the- 
ory (CFT). This CFT ignores short-distance properties 
of the bulk medium, and therefore will not account cor- 
rectly for the term in the entropy proportional to L, but 
it should reproduce correctly the universal constant term. 

To compute the entropy for the case of a disk that con- 
tains an anyon with charge a (far from the boundary), 
we evaluate the partition function Za — tiaS^^^ for the 
associated conformal block of the CFT. Za can be ex- 
pressed as a path integral on a torus of length /3 in the 
Euclidean time direction and length L in the spatial di- 
rection, in the presence of a Wilson loop carrying anyon 
charge a that winds through the interior of the torus in 
the timelike direction. After a modular transformation, 
we have 

Za^Y.^'Zb, (16) 

b 

where Zj, is the partition function for the b block on a 
torus of length L in the Euclidean time direction and 
length /9 in the spatial direction, and S is the modular 
S-matrix of the CFT, which matches the topological S- 
matrix of the anyon model. In the limit L —^ 00, the sum 
is dominated by the trivial block Zi , and we find 



log Za « log (SIZ,) « log 51 + ^{c + c)LIP , (17) 

where c and c are the holomorphic and antiholomorphic 
central charges of the CFT, and S\ — da/T) is a topo- 
logical S-matrix element. Applying the thermodynamic 
identity S = -dF/dT (where F = -TlogZ is the free 
energy), we then find 

Sip) = -^ {T log Z) = aL- log {V/da) . (18) 

Thus when a is the trivial charge and rfa = 1, we recover 
the result of eq. (1) and (2). While this derivation is not 
on so firm a footing as the derivation leading to eq. (15), 
it is more transparent and it generalizes readily to the 
case where the disk contains an anyon. 

We have found an intriguing connection between en- 
tanglement entropy and topological order in two dimen- 
sions. We note that there are close mathematical ties 
between the topological entanglement entropy and the 
(l-l-l)-dimensional boundary entropy discussed in [8], 
and we expect that further insights can be derived from 
studying higher-dimensional analogs of 5topo. We also 
hope that our results can provide guidance for the im- 
portant task of constructing explicit microscopic models 
that realize topological order. 

Results similar to ours have been obtained indepen- 
dently by Levin and Wen [9]. 
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